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EXECUTIVE  SUMMARY 


OBJECTIVE 

The  objective  of  this  work  was  to  derive  the  solution  to  a  static 
ferromagnetic  problem  that  included  both  current-carrying  coils  and  linear 
ferromagnetic  bodies.  The  solution  is  intended  for  comparison  with  the 
solution  to  the  ferromagnetic  problem  obtained  by  various  numerical 
techniques  such  as  the  finite  difference  method,  the  finite  element  method- 
and  the  integral  equation  iterative  solution  method. 

APPROACH 

After  deriving  the  governing  differential  equation  from  Maxwell’s 
equations  for  classical  magnetostatic  field  theory,  the  method  of  separa¬ 
tion  of  variables  was  employed  to  obtain  the  problem  solution. 

RESULTS 

The  magnetic  induction  was  calculated  £or  a  geometry  consisting  of  a 
ferromagnetic  prolate  spheroidal  shell  and  with  a  current-carrying  con¬ 
ductor  both  internal  and  external  to  the  shell.  The  ferromagnetic  body 
was  assumed  to  be  linear  and  homogeneous.  The  reduction  of  the  solutions 
to  that  of  a  prolate  spheroidal  current  band  in  free  space  is  shown  when 
the  permeability  of  the  ferromagnetic  prolate  spheroidal  shell  is  allowed 
to  approach  that  of  free  space  and  when  the  current  in  the  internal  or 
external  prolate  spheroidal  band  is  set  equal  to  zero. 

RECOMMENDATIONS 

It  is  recommended  that  the  derived  solution  be  programmed  on  a 
digital  computer  for  direct  comparison  of  these  results  to  those  obtained 
by  various  numerical  methods.  There  are  plans  to  implement  these  recom¬ 
mendations  during  Fiscal  Years  1980  and  1981. 


ABSTRACT 


The  magnetic  induction  is  calculated  for  the 
configuration  consisting  of  a  ferromagnetic  prolate 
spheroidal  shell  with  internal  and  external  infin¬ 
itesimally  thin  prolate  spheroidal  current  bands. 

The  ferromagnetic  body  is  assumed  to  be  linear  and 
homogeneous.  The  reduction  of  the  solutions  to  that 
of  a  prolate  spheroidal  current  band  in  free  space  is 
shown  when  the  permeability  of  the  ferromagnetic  pro¬ 
late  spheroidal  shell  is  allowed  to  approach  that  of 
free  space  and  when  the  current  in  the  internal  or 
external  prolate  spheroidal  band  is  set  equal  to  zero. 

ADMINISTRATIVE  INFORMATION 

This  ifgrk  was  performed  under  Program  Element  1122 IN,  Project  B005, 
Task  Area  B(^005-SL001,  Work  Unit  1-2704-110.  The  Project  Director  is  Mr. 

W.  J.  Andahazy,  David  W.  Taylor  Naval  Ship  Research  and  Development  Center. 

INTRODUCTION 

As  reported  previously , ^ exact  analytical  solutions  of  Maxwell’s 
equations  using  classical  formulations  have  been  limited  to  body  shapes  and 
inhomogeneities  that  conform  to  a  few  separable  coordinate  systems.  The 
large  computational  and  storage  capabilities  of  modern  computers  permit 
many  electromagnetic  field  problems  to  be  solved  by  using  a  numerical 
solution  to  the  governing  differential  or  integral  equations  under  a 
suitable  choice  of  boundary  conditions.  The  numerical  solutions  of 
Maxwell’s  equations,  when  used  with  a  complete  description  of  the  electric 
and  magnetic  sources  and  the  constitutive  laws  of  the  media,  can  be  used 
to  describe  completely  the  electric  and  magnetic  fields  produced  by  the 
source,  including  nonsyraraetric  geometries,  non symmetric  source  distribu¬ 
tions,  and  spatially  varying  media  parameters. 

The  motivation  for  this  work  arose  out  of  the  need  for  solutions  to 
static  ferromagnetic  problems  that  could  be  used  for  comparison  with  numer¬ 
ical  methods. 


*A  complete  list  of  references  appears  on  page  59. 


PROLATE  SPHEROIDAL  COORDINATE  SYSTEM 
The  prolate  spheroidal  coordinate  system  can  be  formed  by  rotating 
the  two-dimensional  elliptic  coordinate  system,  whose  traces  in  a  plane 
are  confocal  ellipses  and  hyperbolas,  about  the  major  axis  of  the 
ellipses. 

Flammer^  notes  that  it  is  customary  to  make  the  z-axis  the  axis  of 
revolution.  Figure  1  depicts  the  three-dimensional  prolate  spheroidal 
coordinate  system. 


PROLATE  SPHEROIDS,  n  =  CONST 
HYPERBOLOIDS,  6  =  CONST 
MERIDIAN  PLANES,  ip  =  CONST 


e  =  CONST 


if  =  CONST 


Hi 


II  !l 


=  7r/2 


i  =  CONST 


Figure  1  -  Prolate  Spheroidal  Coordinate  System 
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In  this  case,  the  coordinate  surfaces  are:  prolate  spheroids  for  n  =* 
constant;  hyperboloids  of  two  sheets  for  9  =  constant;  meridian  planes 
for  ^  =  constant.  The  prolate  spheroidal  coordinates  shown  in  Figure  1 
are  related  to  rectangular  coordinates  by  the  following  transformation 
equations : 


x  =  a  sinh  r\  sin  0  cos  ip 

y  =  a  sinh  r)  sin  9  sin  ^ 

z  =  a  cosh  n  cos  0  (1) 

where 

0  00 

0  £  IT 

0  2tt 

We  have  denoted  the  interfocal  distance  by  2a  and  the  prolate  spheroidal 
coordinates  by  (r|,0,^). 


BASIC  EQUATIONS 

We  can  start  with  Maxwell’s  equations  for  classical  magnetostatic 
field  problems 


t 


V  x  H  =  J 

V  •  B  =  0 

where 

H  =  Magnetic  field  intensity  (A/m) 

—  2 
B  =  Magnetic  flux  density  (T  or  Wb/m  ) 

—  2 
J  =  Electric  current  density  (A/m  ) . 


(2a)* 

(2b) 


*The  del  operation,  V,  is  defined  with  respect  to  the  rectangular 
coordinate  system  and  is  strictly  valid  in  a  rectangular  coordinate  system 
only.  Very  often,  Vx  and  V*  are  used  as  equivalent  symbols  for  curl  and 
divergence  generally.  This  use  is  followed  in  this  report. 


In  the  general  case  for  ferromagnetic  materials,  B  is  a  nonlinear  function 
H 


B  =  f  (H)  (3) 

where  B  is  not  a  single  valued  function  of  H.  The  function  f (H)  depends  on 
the  magnetic  history  of  the  material;  that  is,  how  the  material  attained 
its  magnetization.  This  is  referred  to  as  hysteresis.  It  is  also  noted 
that  any  property  of  a  ferromagnetic  material  has  meaning  only  if  it  is 
considered  together  with  its  complete  magnetic  history. 

In  certain  practical  engineering  problems,  the  variation  in  the  mag¬ 
netic  intensity  is  small,  and  the  functional  relationship  between  B  and  H 
is  approximately  linear. ^  For  the  linear  case  where  the  material  is 
isotropic,  the  magnetic  induction,  B,  is  related  to  the  field  intensity, 

H,  by  the  relationship 

B  =  y0  (Xnj  +  !)  H  =  Uoyr  H  =  yH  (4) 

=  Magnetic  susceptibility  (dimensionless) 

=  Magnetic  permeability  (H/m) 
y^  =  Relative  permeability  (dimensionless) 

=  Free  space  permeability  (4tt  x  10”^  H/m) . 

This  report  assumes  that  the  ferromagnetic  body  has  isotropic  and 
linear  material  properties.  The  divergenceless  nature  of  the  magnetic 
flux  density  in  conjunction  with  the  fact  that  the  divergence  of  the  curl 
of  any  vector  function  is  zero  allows  the  introduction  of  the  magnetic 
vector  potential  field  A, 

B  =  V  x  A  (5) 

where  A  is  the  magnetostatic  vector  potential  function  in  Wb/m.  The  sub¬ 
stitution  of  Equation  (5)  into  Equation  (2a)  gives  the  fundamental  equa¬ 
tion  of  the  vector  potential  of  the  magnetostatic  field. 


where 


y 

<x«+1> 


(6) 


“V  x  (VxA)  -  (VxA)  x  V  —  =  J 
M  P 

For  homogeneous  materials,  as  assumed  in  this  report,  the  magnetic  perme¬ 
ability  is  spatially  invariant.  Hence, 

vi-0  (7) 

and  Equation  (6)  reduces  to 

V  x  (V  x  A)  *  pJ  (8) 


Using  the  vector  identity 


VxVxA  =  V(V.A)-$A 


(9) 


Equation  (8)  becomes 


V  (V  .  A)  -  $  A  =  yJ 


(10) 


The  magnetic  vector  potential  is  characterized  by  the  important  property 
that  its  divergence  can  be  conveniently  chosen  to  be  zero. 


V  .  A  =  0 


(ID 


Equation  (10)  reduces  to  the  vector’s  Poisson’s  differential  equation. 


$  A  =  -pj  (12) 

This  is  the  governing  equation  for  our  calculations. 

The  general  boundary  conditions  to  be  satisfied  at  the  interfaces  of 
stationary  dissimilar  media  may  be  derived  from  the  limiting  integral 
forms  of  Maxwell’s  equations  and  are  given  by 

"12  '  (i2  '  V  ’  0  “  Bnl  ■  B„2  <l3»> 


5 


V 


(13b) 


n12  X 


(H, 


Hl>  =  Js  or  Ht2 


‘tl 


J 

s 


where  the  subscripts  1  and  2  indicate  the  media  under  consideration  and 
n^2  denotes  the  unit  normal  vector  to  the  interface  and  is  directed  from 
Medium  1  into  Medium  2.  In  the  case  where  the  materials  are  linear  and 
isotropic.  Equations  (13a)  and  (13b)  become 


n12  '  ^M2  H2  “  P1  *V  =  0 


(13c) 


(13d) 


J  is  a  true  surface  current  density  that  may  exist  at  the  interface.  At 
s  _ 

an  interface  where  is  0,  Equations  (13b)  and  (13d)  need  to  be  modified 
accordingly . 


FERROMAGNETIC  PROLATE  SPHEROIDAL  SHELL  WITH  INTERNAL  AND  EXTERNAL , 
INFINITESIMALLY  THIN,  PROLATE  SPHEROIDAL  CURRENT  BANDS 

GENERAL  SOLUTION 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic 
prolate  spheroidal  shell  of  homogeneous  permeability  with  internal  and 
external,  infinitesimally  thin,  prolate  spheroidal  current  bands  of 
constant  current  density  and  ,  respectively.  The  geometry  of  the 
problem  suggests  that  a  prolate  spheroidal  coordinate  system  as  shown  in 
Figure  1  be  used  in  the  problem  solution.  Figure  2,  a  cross  section  of 
the  problem  geometry,  identifies  the  five  regions  of  interest.  The 
boundaries  of  the  prolate  spheroidal  shell  are  determined  by  H  -  and 
n  =  r\  constants.  The  direct  currents  lie  in  the  boundaries  n  =  and 
n  =  r\^9  constants.  Regions  I,  II,  IV,  and  V  have  a  permeability  equal  to 
free  space,  y^,  which  for  convenience  will  be  labelled  y^.  Ampere’s  Law 
states 


V  x  H  -  J 


(14) 


6 


TT* 


n1  ji  vsnrr 


and  since  V  •  B  =  0,  the  induction  B  must  be  the  curl  of  some  vector  field 
A.  The  governing  differential  equation  for  A  when  homogeneous  and  linear 
materials  are  considered,  is  from  Equation  (12) 


0A  =  — uJ 


(15) 


J#) 


CURRENT 

BAND 


CURRENT 

BAND 


Not* 

(  *  eo*h  i) 
v  =  co*  0 


figure  2  -  Cross  Section  of  Ferromagnetic  Spheroidal  Shell 
Surrounding  and  Surrounded  by  Infinitesimally 
Thin  Current  Bands 


The  general  expression  in  prolate  spheroidal  coordinates  for  a 
current  density  is 


J 


Vn  +  Jeee  +  V * 


(16) 


In  the  problem  presented  herein,  the  current  densities  have  only  a  psi  (40 
component  (j^  (6)  e^] ,  which  means  that  the  vector  potential  has  only  a 
psi  component  A(^e^.  The  vector  potential  [  A  =  is  a  function  of  the 


7 


prolate  spheroidal  coordinates  n,  0  [i.e.,  -  A^  (n,9)].  The  constant 

current  densities,  which  lie  on  the  boundaries  between  Regions  I  and  II 
and  Regions  IV  and  V,  can  be  expressed  by  the  functions 


o,  if  0  <  e1  or  0  >  e2 


ve>  V  1E  61 s  9 "  9 2 


where  J^(0)  is  equal  to  a  constant  along  v\  *  for  0^  £  6  £  0^  and 


0,  if  0  <  0^  or  0  >  02 


l  Ve)  V  lf  e!s  0  s  9 2 

where  J^(0)  is  equal  to  a  constant  J2  along  n  *  for  0*^  £  0  S  6^. 
Therefore,  Equation  (15)  has  only  an  azimuthal  or  psi  component  and  can 
be  expressed  as 

O  A  =  (n,0)  =  0  in  Regions  I  through  V  (1 

When  the  vector  Laplacian  is  expanded  in  prolate  spheroidal  coor- 

dinates,  Equation  (19)  can  be  expressed  (see  Appendix  A,  reference  2) 


r_i 

I  sinh 


3  (sinh  r|  A^ 

n*  9n 


ih-iU 

J  30  [sin 


9  (sin  0  A^ 
0  30 


=  0  (20) 


Applying  the  method  of  separation  of  variables,  let  us  assume  that  A  can 
be  expressed  as  the  product  of  two  functions 

A,  =  H  (cosh  n)  G  (cos  0)  (21 


where  H  (cosh  n)  is  a  function  of  cosh  n  only,  and  G  (cos  9)  is  a  function 
of  cos  0  only. 

Substituting  this  form  of  the  vector  potential  A  into  Equation  (20) , 
we  have  after  separation  of  variables 


^-4r  +  coth  n  ^  -  (v  (P+1)  + - “2 

dn2  n  \  sinn  n 


H  =  0 


(22a) 


+  cot  9  ^  +  ^p  (p+1) 


sin2  9 


G  =  0 


(22b) 


where  the  separation  constant  is  p  (p+1)  and  p  is  an  integer  from  one 
to  infinity.  It  is  well  known  that  differential  equations  of  the  form 


2  „  ^  /  2 
H  H  .  dH  /  m 

+  coth  n  -T- —  p  (p+D  +  o 

dn2  dn  \  sinn  n 


H"  =  0 


(23a) 


have  the  general  solution  of  the  form 


H'  =  ClPm  (cosh  n)  +  C0Qm  (cosh  n) 

1  p  2  p 


(23b) 


where  and  are  constants  and  it  is  known  that  a  differential 
equation  of  the  form 


d2G'  ,  ^  n  dC' .  /  ,  m2 

~~2~  +  cot  9  dO  +  P  (P+1)  ‘  “TT 

dO  \  sin  0 


G>  =  0 


(24a) 


has  the  general  solution  of  the  type 


G'  =  C_Pm  (cos  0)  +  C,Qm  (cos  9) 
3  p  4  p 


(24b) 


where  and  C,  are  constants.  Pm  and  Qm  are  the  associated  Legendre 
3  4  p  p 

functions  of  the  first  and  second  kind,  respectively. 


Comparison  of  Equations  (22),  (23),  and  (24)  shows  that,  in  Equations 
2 

(23)  and  (24),  m  is  equal  to  1.  This  requires  that  m  always  equals  unity. 
The  solutions  of  Equations  (22a)  and  (22b)  are  expressed  as 

H  (cosh  n)  =  Apl  (cosh  n)  +  BQ^  (cosh  n)  (25a) 

P  P 

G  (cos  0)  =  A'pJ  (cos  0)  +  B'qJ  (cos  0)  (25b) 

The  general  solution  of  Equation  (20)  may  be  formed  from  the  product 
of  solutions  in  Equations  (25a)  and  (25b) ,  which  yield 


=  H  ,(cosh  n)  G  (cos  9) 


(cosh  n)  G  (cos  0) 
P 


(26) 


A.  =  I  AP1  (cosh  n)  +  BqI  (cosh  X]) 

V  i  I  P  P 

p=l  L 

x  j^A^Pp  (cos  0)  +  (cos  0)  J 


(27) 


pm  and  Q  are  the  associated  Legendre  functions  of  the  first  and  second 
P  P 

kind,  respectively. 

For  the  prolate  spheroidal  system,  the  associated  Legendre  functions 
of  the  second  kind  are  infinite  at  cos  0  =  ±1,  and  as  such  cannot  be 
included  in  a  general  solution  for  a  given  region  which  includes  0  =  0  or 
0  =  it.  Therefore,  in  our  case,  the  constant  is  set  equal  to  zero. 
Equation  (27)  reduces  to 


OO 

■  r  vi 

p=l  L 


(cosh  n)  +  'i  (cosh  n)  P  (cos  0) 


1  p1 
P 


10 


(28) 


When  the  substitu- 


where  and  are  constants  (K^  =  AA' ,  =  BA  ) . 

tions  £  =  cosh  n  and  V  =  cos  0  are  made  in  Equation  (28),  A^  can  be 
expressed  as 


\  ■  E  [Vj®  +  v>J®  1  ei 

P=1  L  J 


(V) 


(29) 


This  is  the  general  form  of  the  psi  component  of  the  vector  potential 
that  will  be  used  to  determine  the  potentials  A^  in  each  direction. 

BOUNDARY  CONDITIONS 

The  form  of  the  component  of  the  vector  potential  A^  in  Regions  I 
through  V  is  determined  from  Equation  (29).  These  magnetostatic  vector 
potentials  in  Regions  I  through  V  are: 


n 


>ii 


■  t  [vH 

P-l  L  J 


pj(v) 


■  E  IVp®  +  cPQJ®  1  pp1(v) 

p=l  L  J 


>111 


■  £[vJ®  +  v$®H 

p=i  J 


(V) 


(30) 


>IV 


oo 

-  E 


P=1 


F  P  (O  +  GQ 
P  P  P  P 


;®  ] 


(V) 


>v 


■  Efvi®l  plpM 

p^l  L 


Because  the  potential  must  be  finite  in  each  of  the  Regions  I 
through  IV,  and  approach  zero  as  £  **  00  in  Region  V,  the  following 
constants  were  set  equal  to  zero: 

11 


a.  For  A  ,  the  constant  associated  with  Q^(£)  P^(V)  was  set  equal 

uu  pp 


to  zero  because 


Qp(£)  -+  00  at  £  =  1  (z  axis  between  ±a) 

b.  For  A  ,  the  constant  associated  with  P^(£)  P^(V)  was  set  equal 
ipv  p  p 

to  zero  because 


Pp(0  00  as  £ 


(we  note  Q  (£)  0  as  £  °°)  • 

P 

A,B,C,D,E,F,G,  and  H  are  constants  to  be  determined 
PPPPPPP  P 

from  the  boundary  conditions.  At  each  interface,  the  basic  laws  of  mag¬ 
netostatics  (Equations  (2a)  and  (2b))  reduce  to  boundary  conditions  on 
B  and  H  that  can  be  used  to  evaluate  these  eight  constants.  The  normal 
component  of  B  across  each  boundary  must  be  continuous;  i.e.,  (B2  ~  B^) 
n12  =  0  where  the  quantity  n^  is  the  unit  outward  normal  to  the  surface. 
This  provides  the  following  boundary  conditions  which  must  be  satisfied 
by  the  solutions  given  in  Equation  (30)  for  each  region 


Bni  =  Vi  n  =  ni 


(31a) 


Vi  =  Vn  n  =  n2 


(31b) 


Vn  =  Vv  n  =  n3 


(31c) 


Vv  "  Bnv  n  =  n4 


(31d) 


The  eta,  or  normal,  component  of  the  magnetic  field  is  expressed 
in  terms  of  the  vector  potential  as 


3  (e»A, ) 


B  =  (V  x  A ,  )  = - - 57T 

H  '.p  r\  e2e3  30 


V  55  \  }  <32> 
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where 


B  =  V  x  A  = 


2  2 

a(sinh  n  +  sin  0)(sinh  H  sin  6) 


_  O  <y  1  _  n  ft  i  _ 

(sinh  n  +  sin  0 )  2  eg  (sinh  f|  +  sin  6)  2  e^(sinh  n  sin  6)j 


a_ 

3n 


8 

90 


9_ 

9i^ 


sinh  n  sin  0 


and 


2  o  k  2  2k 

£  =  cosh  n,  ex  =  e2  =  a(sin^  H  +  sin  0)^  =  a{£  -V  )  ,  V  =  cos  9 
■  a(sinh  n  sin  0) 

However,  since  the  vector  potentials  in  each  region  are  functions  of 
P1(v),  we  can  simplify  Equation  (32)  to  constraints  on  at  the 
interfaces : 


A ,  T 

\ll 

at 

li 

JZf 

(33a) 

=  AlpIII 

at 

n  =  n2 

(33b) 

NjIII 

> 

M 

11 

at 

n  =  n3 

(33c) 

A<piv 

~  Aipv 

at 

n  =  n4 

(33d) 

The  second  set  of  boundary  conditions  states  that  the  theta,  or 
tangential,  component  of  H  across  each  boundary  must  satisfy  the  rela¬ 
tionship 


12 


(H2  -  V  =  Js 


(34) 


i 


13 


where  Jg  (which  equals  J^(0))  is  the  real  surface  current  density  in  the 
limit  of  vanishing  width  between  the  two  regions.  Using  the  relationship 
B  =  pH,  Equation  (34)  can  be  expressed  as 


!s!2  .  hi  .  .  ... 

v2  Pj  * 


The  current  must  be  expanded  in  a  series  of  associated  Legendre 

functions  P^(v)  as  in  Reference  2.  The  form  of  the  current  is 
P 

00 

J  £  Vp  (cos  6) 

J,  (Q)  =  — p-x2  - 2 — r 

a(sinh  r\  +  sin  0)^ 


where  using  t,  =  cosh  U,  V  =  cos  0,  can  be  shown  to  be 


=  - (2p+l)  a 
p  2p(p+l) 


(C  -V  ) 2  P  (V)dv 


For  the  two  current  bands  of  interest,  we  have 


V0) 


ji  £  yfo> 

a(^~V2)h 


(38a) 


where 


=  - (2p+l)  a 

p  2p (p+1) 


/ 2  * 


(v)dv 


(38b) 


and 


The  general  expressions  for  the  potentials  in  each  region  (Equation 
(30))  are  then  substituted  into  the  boundary  conditions  (Equations  (33) 

and  (41))  and  solved  for  the  eight  constants  (A,B,C,D,E,F,G, 

^  P  P  P  P  p’  p’  p’ 

and  Hp) .  Since  there  are  eight  equations  with  eight  unknowns,  the  poten¬ 
tial  in  each  region  can  be  specified.  The  eight  boundary  value  equations 
are  presented  below.  The  index,  p,  in  the  summation  sign  has  both  even 
and  odd  value  and  it  takes  on  values  from  1  to  The  eight  expressions 
for  the  boundary  conditions  are: 

vJ<5i)pp1(v)  -[vJ(5i)  +  %ql/v  ]  p>>  <«»> 


k  [‘S2-1)1*  |»ppJ<5)  +  cpQPCE)|p\v>] 

C-«l 


U1  I  a^-v2)*5 


+ 


a(^-v  ) 


d_ 

U 


(yj«>)pp1<v)] 


(8) 


(42b) 


Vp«2>  +  Cp<1P<52>]  PJ(V)  '[VJ'V  +  V^V]  rl 


(V) 


(42c) 
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V 


1  \  1  9 


VattJ-v2)*  H 


(C2-!)**!  HpQp(^)  JPp (V) 


1  \  1  9 


"i  LkW)1*  8C 


(C2-!)15  FppJ(C)  +GpqJ«)  ]pJ(v) 


J  ,(6) 


(42h) 


By  making  the  following  substitutions 


pp«>  -!i[«2-1>'5  pJ(5)] 


and  performing  algebraic  manipulation,  the  eight  boundary  conditions  can 
be  simplified  to: 


A  P1(C1)  -  B  phs,)  +  c  <£«,) 
PP  I  PPI  PP1 


(43a) 


—  ^  [b  PA(£  )  +  C  QA(?i  )  ]  +  / ■—  ^  f A  pA(C1)l 

ji y L p p  1  p p  1 J  \yi/Lp p  'J 


Jpl(0)  a(^-V2)h 

Phv) 

P 


(43b) 


The  solution  of  these  eight  simultaneous  equations  to  obtain  the  constants 
gives : 


A  =  B  +  C  -f — —  (44a) 

P  p  p  Pi(5i) 


c 

p 


W1-Jpl(e)  a(^-v2) 


pj(v) 


PA(CO 

p  1 


Qp«l> 

L  Pp^l^ 


(44b) 


F 

P 


fe  _e 

pJ<53)  P  pJ(C3) 


+  D 


(44c) 


E 

P 


■  Dp  [Q]  -  %  DO 


(44d) 


where 


[o]  ■ 


«J«3> 


*J«3> 


W  pj({3> 


[R]  - 


<<t3> 

P^3> 


£«3> 


pJ«3> 


LPJ«3> 
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/yl\  Qp(^3> 

\y2/PA(C.) 
P  ^ 


Since  the  eight  coefficients  can  be  determined  for  a  specified  problem 
from  Equation  (44),  the  potentials  Aj. ,  A^,  AIV’  anc*  in  ReSions  * 

through  V  can  be  completely  determined.  The  normal  (B^)  and  tangential 
(Bq)  components  of  the  magnetic  induction  in  each  region,  I  through  V,  can 
be  determined  by  using  Equations  (32)  and  (40). 


APPENDIX  A 


DERIVATION  OF  COEFFICIENTS  OF  THE  VECTOR  POTENTIAL  FOR  A 
FERROMAGNETIC  PROLATE  SPHEROIDAL  SHELL  WITH  INTERNAL  AND 
EXTERNAL,  INFINITESIMALLY  THIN,  PROLATE 
SPHEROIDAL  CURRENT  BANDS 


In  this  appendix,  the  coefficients  are  derived  for  the  vector  poten¬ 
tial  in  Regions  I  through  V  for  a  ferromagnetic  spheroidal  shell  with 
internal  and  external,  infinitely  thin,  spheroidal  current  bands.  For  a 
detailed  discussion  of  the  ferromagnetic  problem,  see  the  body  of  this 

report.  The  components  of  the  magnetic  vector  potential  in  each  region 
is  given  by 


A,  T 


ZJa  p\t) 

p=i  p  p 


PJ(V) 


(A. la) 


A(|/II 


2  [b  +  C  qj(5)l  ?hv) 

P=1  P  P  P  P  J  P 


(A. lb) 


Vi i  '  £  (Vi®  +  epqp(5)]  pp(v) 

oo 

Vv  ■  P5  [vJ® +  gpqp1(5)]pJiv) 


CA.lc) 


(A. Id) 


E 

P=1 


vJ(5) 


pj<v) 


(A.le) 


The  coefficients  in  Equations  (A. la)  to  (A.le)  are  obtained  by  substi¬ 
tuting  these  equations  into  Equations  (33)  and  (41). 
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Vp({l>  +  r> 


(=r)  irt?  £  [(52-1)! 


B  p‘(0  +CQ1 
p  p  p  p 


H 


_1 _ 9 

M  H 


k2-1)|5  (yj«>)  p>> 


'p1<£2>  +  CpQJ<{2)]  ’’p1<V)  ‘  [DpPp1<52>  +  Vj«2>]  Pp1(v> 


at^-v2)*  S« 


K2-1)'i  (dpJ(C)  EpQJ(C)j  pj(v) 


a(d-v2)li 
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vJ<o  +  cpqJ<E)  PJ<U) 


n«]M  *'**?m»F 


•  t . 


(fj  si  [<s2-1)!i  +  v$<0)  pp 


(v) 


f=.f 


'3 


2  2.k  ac 


a(f,:-v  ) 


(^-D15  (yj(0  +  *pqJ<«)  pJ(v). 


'3 


(A. 2f ) 


[fppp<^>  + 


PJ(V)  ’  Hp«J(V  PJ(V) 


(A.2g) 


('  *1 )  a(f4-v2)’5  3« 


a2-Dh  (vi<£))pJ<w) 


+ 


,_2  2. 
a(£,-v  ) 


_a_ 


"<e2-i>‘s  (> 


F  P^(0 
P  P 


G  Q1 (C) 
P  P 


PL(v) 

P 


Jp2(0> 


(A. 2h) 
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V 


Define 


Pp(0  =  g|  UC2-!)4  Pj(01 


QA(£)  =  [(C2-!)1*  Q*(0] 

p  ^  p 


(A. 3a) 


(A. 3b) 


After  appropriate  substitution  of  Equations  (A. la)  to  (A.le)  into  Equations  ’■ 

(33)  and  (41)  and  using  Equations  (A. 3a)  and  (A. 3b),  the  following  £ 

boundary  value  equations  are  obtained.  [.: 


A  P1(51)  =  B  P1(51)  +  c  Q1(C1) 
ppl  ppl  ppl 


p  p  1  p  p 


(A. 4a) 


i)  (bppp«i>  +  cpqp«i)) + (i)  (v^v 


Jpi(9)  aC^-v2)^ 

PL(v) 

P 


(A. 4b) 


VJ«2>  +  CpQi«2>  -  vi(«2>  +  Ep<<52> 


(A. 4c) 


i)  [DPPP(E2> 


e  <£(£,)]  = 

P  P  2 


—  |  [B  PA(C„)  +  c  QA(CJ] 
,Pj  PP  2'  p  P  2 


(A.4d) 
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iftf 


(A. 4e) 


3  P1(C.)  +  E  Q1(C.)  =  F  pj(5.)  +  G  qV  .) 
p  p  3  p^p  J3  p  P  3  P  P  3 


+  EPQPa3>i  - 


—  )  [F  pA(£_)  +  «<£(£,)  J  (A .40 


M  p p  3  pp  3 


F  P1^,)  +  G  Q1  (C/ )  =  H  QJ  (C,) 
p  p  4  p  p  4  p  p  4 


(A.  4r.) 


^)ihpqp(V>  +(^)[Vp'V  +  0pqX» 


Jp2(0)  a (g-v  )'2 

pl(v) 

P 


(A. Ah; 


These  algebraic  equations  provide  eight  simultaneous  equations  with 

eight  unknowns;  and  they  can  be  solved  for  the  coefficients  A  ,  B  ,  C  , 

~  P  P  P 

Dp,  Ep,  Fp,  0^,  and  H  .  Solving  Equation  (A. 4a)  for  A^  gives 


A 

P 


qp1(5i> 

B  4*  C  -f - 

p  P 


p  w 


The  solution  of  A 

P 


from  Equation  (A. 4b)  is 


A 

P 


B  +  C 
P  P 


pp<V 


2  2  k 

J  1 (0)  a(f"-v fc) 2 

-El _ A _ 

phv)  pA(cp 

p  pi 


(A.  b) 


(A.  6) 
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JUl 


\ 

\ 

■i 


Equating  Equations  (A. 5)  and  (A. 6)  and  solving  for  C  gives 

P 


Uj,  J  x(0)  aCC^-v2)*5 


p  pJ(v)  P^(C,) 


pJ«i>  pp«i> 


(A.  7) 


C  =  J\ 
P  Pi 


(A. 8) 


Solving  Equation  (A.4e)  and  (A.4f)  for  gives 


Q*(5J  (£(£,) 

F  =  E  -? — —  -  G  -2 — —  +  D 

p  p  P*(£  )  P  PX(CJ  P 
p  3  p  3 


(A.  9) 


r  (H)  *  VjS«3>l  c  <«3> 

»  '  \v2/  p‘(53)  p  pJ(53) 


Equating  Equations  (A. 9)  and  (A, 10)  and  solving  for  E^  gives 


where 


tQi = ^  -1 


Ep  -  “p  to)  -  Gp  IR1 


_  /M 

pJ(C3)  W  p‘«3> 


(A. 10) 


(A. 11) 


(A. 12) 


LW  WJ/  LPp(C3> 


V  2  /  Pp«3)J 


E 

P 


Ly2 


-1 


.  o  +  G.  Q^> 


p  p„A(5,)  P  Pi<5,> 


p-*3'  *  p  vv*3 


/m  QXi 


PJ«3> 


‘2/  P“«3) 


(A. 14 


Solving  Equation  (A. 4c)  for  Bp  and  substituting  Equations  (A. 8)  and  (A. 11) 
gives 


B  =  J1?  +  D  [S]  -  C  [Z] 


P  Pi  P 


(A.  15 


where 


,11 


j  ;  =  -  j. 


i 


Pi 


Pl  PJ«2> 


(A. 16 


[S]  = 


Qp«2> 
l  +  [Q]  f  — - 


(A  .17 


The  solution  of  Equation  (A.4d)  for  and  using  Equations  (A.8)t  (A. 11), 
and  (A. 15)  gives 


(srKi  qp«2>  p-k’> 


Pi  P  2 


D  = 


('  S^)  151  Pp«2>  *(4)Pp«2>  +(^)W  <«2> 


+  G 


(^)1Zi  PP«2>  <«2> 


[K)'S1  PP<?2>  +(^)PpK2>  +(^)'QI  Qp<52>"" 


(A. 19) 


where 


I1?1 

Pi 


[U]  = 


D  =  J1*1  +  [U]  G 
P  Pi  P 


j\  qa(cj  +  4V1;  pa(^9) 

,  Mi  J  Pi  P  2  Mx )  pl  P  2 


[  ^)tSl  Pp«2>  +(^)Pp<52>  *(V2)M 


K) 


«7)  [Z!  PPC52>  + 


H)  [R1  °P«2> 


.('  *1 ) 


ttM  [S]  PpU2)  + 


KK'v  *K) 


IQ]  Qp(C2) 


(A. 20) 


(A. 21) 


(A. 22) 


Now,  solving  Equation  (A.4g)  for  F  gives 

P 


F  *  H 
P  P 


[A] 


(A. 23) 


where 


[A] 


<>x> 

pJ<e4> 


(A. 24) 


The  solution  of  Equation  (A.4h)  for  F  gives 


F 

P 


Gp  [Cl  +  Hp  [C|  +  jJ2 


(A. 25) 


where 


[Cl 


qX> 

pp(v 


(A. 26) 


Jp2(9)  a(^-V2)4 

pp«4)  pj>> 


(A. 27) 


Equating  Equations  (A. 23)  and  (A. 25)  and  solving  for  H  ,  gives 


where 


H 

P 


G 

P 


+  J 


II 

P2 


(A. 28) 


(A. 29) 
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Now,  solving  Equation  (A.4e)  for  using 


[H]  = 


Pp«3> 


(A. 30) 


gives 


F  =  E  [H]  -  G  [II]  +  D 
P  P  P  P 


(A. 31) 


Now,  solving  Equations  (A. 23)  and  (A. 31)  for  using  Equations  (A. 11), 
(A. 20),  and  (A. 28)  gives 


-j”1  [Q]  [H]  -  j"1  +  [A] 

Jp  IU]  [Q]  [H]  -  [R]  [H]  -  [H]  +  [U] 


(A. 32) 


where 
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«i«3>  (vl\  Qp«3> 


pJ(53)  \p2/  p*(53) 


(A. 33) 


[R]  = 


0p^3> 
Pp(?3>  PJ«3> 


«X>  /pl\  Qp«3> 


pJ(53)  '  2 /  PpCS3) 


(A. 34) 


[S]  =  1+  [q](q;«2)/pp 


>4 


[Z]  = 


[l»l  / PJ(52>] 


(A. 35) 


(A. 36) 


APPENDIX  B 


DETERMINATION  OF  THE  MAGNETIC  VECTOR  POTENTIAL  FOR  AN  INFINITESIMALLY 
THIN  PROLATE  SPHEROIDAL  CURRENT  BAND 


In  this  appendix,  the  potentials  A^  in  the  inner  region  and  A^^  in 
the  outer  region  are  derived  for  the  infinitesimally  thin  current  band  in 
a  homogeneous  medium  of  permeability  y^  (Figure  B.l). 


Figure  B.l  -  Infinitesimally  Thin  Spheroidal  Current  Band 


The  potential  in  the  inner  region  and  the  outer  region  A^j.  of 
the  infinitesimally  thin  current  band  problem  are  solutions  to  the  vector 
Laplace's  equation  0  A  =  0.  These  solutions  can  be  expressed  as: 


rv . 


£  [A  pj(5>]  pj(v) 

P=1 


(B.la) 


>11 


£  [F  qj(5)l  pJ(v) 

p=l 


(B. lb) 


The  coefficients  A  and  F  are  determined  from  the  boundary^ conditions  of 
P  P 

the  problem.  After  algebraic  manipulation  such  as  with  Equations  (32)  and 
(AO)  in  the  text,  the  boundary  conditions  for  the  normal  component  of  B 
and  the  tangential  component  of  H  become: 

Aj  «  A  at  n1  =  (B.2a) 


-  ye) 


(B . 2b) 


Using  the  relationship  £  -  coshr)  and  v  =  cos0  and  substituting  the  expres¬ 
sions  for  A^j.  and  (Equations  (B.la)  and  (B.lb))  into  the  boundary 
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value  equations  (Equations  (B.2a)  and  B.2b))  provides  us  with  the  following 
algebraic  equations  for  the  coefficients 


AP^C,)  =  F  QV  ) 

p  p  1  p  p  1 


(B.3a) 


i  A  .  A  JO)  aCSfV)' 

-  ^  'W^1 +  ^  'Vp'V  ■  — 


(B . 3b) 


where  the  following  substitutions  were  used 


P*({>  -  si  US2-!)1*  pJ(£)J 


(B.4) 


Qp(C)  =  af  [(C2-!)*5  QpCC)  J 


(B.5) 


These  equations  are  solved  for  and 


F  by  simple  algebraic  manipulations 
P 


A 

P 


?  qp(€i) 

PpX> 


(B . 6a) 


F 

P 


«foi>  <«!> 
pJ«i> 


(B.6b) 


The  potential  A.j  and  A^j  are  determined  by  substituting  the  expression 

for  A  and  F  into  Equations  (B.la)  and  (B.lb). 

P  P 
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n 


& 


Jp(9)  a(g2-v2)>i 

»S«i>  pJ<v> 


PJ«1>  pp«l> 


x  ?ho  p hv> 

p  p 


(B.7a) 


(JjII 


E 

p=l 


y.  J  (0)  acd-v2)15 
1  P _ I _ 

Pp^l)  P£^ 


<«i>  _  ^«i> 


LpJ«l>  pp«l> 


qj(5)  pJ(v) 


(B. 7b) 


where 


JG  P  (v) 

JP(e>  =  -Vn 

P  a(^-v2)'5 


(B.8) 


APPENDIX  C 


REDUCTION  OF  THE  MAGNETIC  VECTOR  POTENTIAL  FOR  A  FERROMAGNETIC 
PROLATE  SPHEROIDAL  SHELL,  SURROUNDING  AND  SURROUNDED  BY 
INFINITESIMALLY  THIN  CURRENT  BANDS,  TO  THAT  OF  A  THIN 
COIL  IN  FREE  SPACE  WHEN  THE  INTERNAL  CURRENT  DENSITY 
[J  (6)]  IS  ZERO  AND  IN  THE  LIMIT  y  EQUALS  y 


In  this  appendix,  the  coefficients  A,B,C,D,E,F,G,  and  H  , 

ppppppp  p 

for  the  potentials  are  calculated  for  the  system  consisting  of  a  ferromag¬ 
netic  shell  of  permeability  y^  with  internal  and  external  current  bands  in 
a  homogeneous  medium  with  permeability  y^  when  the  internal  current  density 
[J^(0)]  is  set  equal  to  zero  and  the  limit  is  taken  as  y^  =  y^ .  These 
coefficients  are  utilized  in  Equation  (30)  in  the  body  of  this  report. 

The  variables  are  defined  in  Figure  2  of  the  text.  When  J^(0)  is  set 
equal  to  zero  and  y^  is  set  equal  to  y^,  the  problem  reduces  to  that  of 
finding  the  potential  in  the  two  regions  of  a  simple  current  band  (Figure 
B.l  in  Appendix  B) ,  because  the  ferromagnetic  shell  will  now  have  a  perme¬ 
ability  y^  equal  to  that  of  the  homogeneous  medium  with  permeability  y^. 

In  this  limit,  the  coefficients  should  assume  the  following  form 

A  =  B  =  D  =  F  (C.la) 

P  P  P  P 


C  =  E  =  G  =0  (C.lb) 

P  P  P 

and  H  and  F  should  reduce  to  the  coefficients  for  the  potentials  in  the 
P  P 

two  regions  for  the  prolate  spheroidal  band  problem  (see  Appendix  B) .  It 

the  coefficients  assume  this  mathematical  form,  it  will  prove  that  the 

mathematical  form  of  the  coefficients  for  the  prolate  spheroidal  shell 

with  the  external  coil  are  mathematically  correct. 

The  mathematical  solution  for  G  in  terms  of  known  quantities  was 

P 

derived  in  Appendix  A  and  was  reported  in  the  text  (Equation  (44h)). 

M1  [Q]  [H]  -  M  +  J*2  [A] 

r  =  _ El - - — - — — P-i - kL -  (C  2a) 

P  [u]  [Q]  [H]  -  [R]  [H]  -  [H]  +  [U]  V  ' 
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limit  [U]  lQ]  [HJ 


!  (*V  ^(5,)  | 


w  w 


p;<9 

a,, , 


V%>  Qp(C2)  t  Qg(S2) 

PP1(«3)  <(«3>  PJ«2>  Pp<S2> 


|/q*<£,)  QP<5,)\  //p5i>  <£({,> 

limit  [ R 1  (HJ  -  |(-Pi---P3-)/(-p3 - p1- 

“2^1  IV, <?3>  W//  \Pp'53)  't'P 


_FL_A_ 

p;<5l) 


therefore 


limit  [ U ]  -  oo 
U  1 


1*  i 


Now,  because 


F  =  H  [A]  -  G  [A] 
P  P  P 


where 


F 

P 


H 

P 


[A] 


P2=Pl 


[A] 


pJ<v 


Also 


H 

P 


G 

P 


+  J 


II 

P2 


Therefore 


H 

P 


U2-Vi 


Ml  Jp2<6)  a(t*- 

PX>  P>> 

'&V  <<5 


ll2=Ul 


Hi  Jp2 (6)  a 

pA(s/)  phv) 

p  4  p 


’<>X> 
PJ«4> 


“*  / 
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' .  ■»  * »ji» rt+mammm 


A 


- 


and,  using  Equations  (C.10)  and  (C.22),  we  have 


=  I) 

y1  Jp2(0)  a ^l~v2)h  /qJ(c4)V 

/( 

p 

=yi  'V 

=U1 

.  pp«4)  pp^»  vXv. 

/[ 

Q  (£/) 

pJ(C, ) 

P  4 


qV4) 

PA(v  ) 
P  4 


A  derived  expression  for  B  was 

P 


(C.2A) 


B  =  J11  +  D  [S]  -  G  [Z] 
P  Pi  P  p 


(C. 25) 


d,  using  Equations  (C.3),  (C.6) ,  and  (C.10),  we  have 


B  I  =  D 
p  i  P 

u2=yl  p2=Pl 


(C. 26) 


Thus  we  have 


A  -  B 

-  D 

=  F 

P  P2=ul  P  y2 

M2=yi 

p  *‘2 

"i 

with 

1 

Jp2(0)  aKl-vZ)h  , 

(M] 

/ 

Qp(44) 

p  * 

u2=p, 

pX>  ^  ' 

/ 

Phf,.) 

P  A 

, ) 

P  *  J 

(C. 27) 


y2=y] 


Q  (4,  ) 

JP _ 

.pV,,) 

P  A 


iC.28) 


A7 


and 


C  =  E  =  G  =0 
P  P  P 


(C. 29) 


and 


yl  JP2(0)  a^4"v2>4 

pX>  pJ<'-> 

Pp«4>  Pp«4> 

(C. 30) 


Now #  Equation  (30)  reduces  to 


A\pn  A\pIII  A\piv 


P-1 


.A 

p  ! 

V 

V 

pV)  P*(v)  (C . 31) 
p  p 


oo 

**»■  Jf 


H 

P 

"V 

«J<»  p>> 


(C . 32) 


comparing  Equations  (C.31)  and  (C.32)  to  Equations  (B.la)  and  (B.lb), 
which  are  repeated  for  convenience  as  and  using  a  primed  notation 


n 


-  X  A„  p?(o 

P=i  p  p 


pj<v) 


(C . 33a) 


Fn  <£<5)  PJ(V) 

P  P  P 


where 


A' 

P 


F' 

P 


■fry 

pp«1> 


(C. 33b) 


(C.34) 
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Thus 


Therefore , 


F  = 
P 


(C. 35) 


A'  =  F  I  or  A'  -  A  I  (see  Equation  (C.27) 
P  P  DP1 

VW1  U2=M1 


With  q  -  c4,  JpO)  =  Jp2(0) 


(C. 36) 


¥'  =  H 
P  P 


U2=yl 


with  q  =  C4,  Jp(6)  =  Jp2(0) 


(C. 37) 


in  the  limit 


and  with  J  . (0)  =  0 
Pi 


\ll  \lll  \lV  \l 


(C . 38) 


(C . 39) 
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APPENDIX  D 


REDUCTION  OF  THE  MAGNETIC  VECTOR  POTENTIAL  FOR  A  FERROMAGNETIC 
PROLATE  SPHEROIDAL  SHELL,  SURROUNDING  AND  SURROUNDED  BY 
INFINITESIMALLY  THIN  CURRENT  BANDS,  TO  THAT  OF  A  THIN 
COIL  IN  FREE  SPACE  WHEN  THE  EXTERNAL  CURRENT  DENSITY 
(J  (0)]  IS  ZERO  AND  IN  LIMIT  p,;  EQUALS  p 


In  this  appendix,  the  coefficients  A  ,  B  ,  C  ,  I)  ,  E  ,  F  ,  G  ,  and  H 

p  p  p  p  p  p  p  p 

for  the  potentials  are  calculated  for  the  system  consisting  of  a  ferromag¬ 
netic  shell  of  permeability  with  internal  and  external  current  bands  in 
a  homogeneous  medium  with  permeability  p^  when  the  external  current  density 
[ J 2 O ) ]  is  set  equal  to  zero  and  the  limit  is  taken  as  p ^  *  p^.  These 
coefficients  are  utilized  in  Equation  (30)  in  the  text.  The  variables 
are  defined  in  Figure  2  of  the  text.  When  ^(B)  is  set  equal  to  zero  and 
pn  is  set  equal  to  p^,  the  problem  reduces  to  that  of  finding  the  potential 
in  the  two  regions  of  a  simple  current  band  (Figure  B.l  in  Appendix  B) , 
since  the  ferromagnetic  shell  will  now  have  a  permeability  p^  equal  to  that 
of  the  homogeneous  medium  with  permeability  p^. 

In  this  limit,  the  coefficients  should  assume  the  following  form 


B  «  D  =  F  =  0  ( l) .  1  a ) 

P  P  P 


C  =  E  -  G  =  H  (D. lb) 

P  P  P  P 


and  where  A  and  C  should  reduce  to  the  coefficients  for  the  potentials  in 
P  P 

the  two  regions  for  the  prolate  spheroidal  band  problem  (see  Appendix  B) . 

If  the  coefficients  assume  this  mathematical  form,  it  will  prove  that  the 

mathematical  form  of  the  coefficients  for  the  prolate  spheroidal  shell 

with  the  internal  coil  are  mathematically  correct. 

The  mathematical  solution  for  G  in  terms  of  known  quantities  was 

P 

derived  in  Appendix  A  and  was  reported  in  the  text  (Equation  (44h)). 


c; 

p 


-j 


ill 


.ill 


.ii 


'Pi  IQl  l«]  -  V  +  V  [A' 
[UJ  fQ)  [H]  -  [R]  [H]  -  [Hi  +  [ill 
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U).2a) 


where 


(mT)  JP1  FP(^> 


■  V 


[A]  -  [C] 


[U]  = 


’(-  V  ^(?2>  +  (sj)  'S’  ♦  (V  <«2: 


■fi 


“i  \  / [Qp<e3>  (vj,\i^i 

U2  )  P1^,)  \g2/p\Co) 


Z)  -  [l«l  <iJ<S2)y/pJ<52>] 


(])•*-'  j) 


Uj  j  (0)  a^-V2* 


P  (L)  P  (v) 
p  4  p 


(I).  2k) 


The  coefficient  Cp  will  now  be  evaluated  when  J2(0)  =  0  and  the  limit  is 
taken  with  =  When  =  ® 


j1;  -  o 

p2 


(D.3) 


and  G  reduces  to 
P 


-jpT  [q1  [h]  -  il1 

[U]  [Q]  [H]  -  [R]  [Hi  -  [Hi  +  [111 


(D.4) 


divide  numerator  and  denominator  by  [U]  and  when  U2  '♦  Mj 


limit  pi 

VM1  [Hi 


jh  -  jpi 


<^2> 


<ca,)  p"(c?) 


pno 

p  2 


W  ±rr  -  w 


p 


(II.  5) 


limit  [R]  = 
VJ  2’>VJ  1 


(D.6) 


limit  [Q]  =  0 

Vyi 


(D.7) 


Thus 


Vyl 


limit  [U]  = 

00 

Vyl 

[Q] 

[H] 

limit  — 

-  =  0 

Vyi 

[U] 

Jpl(0)  a(^-v2)% 

/ 

«x>  <«i> 

p‘(v)  pJtSj) 

/ 

pp(ei'  Pp«l> 

(D.8) 


(D.9) 


=  Jpl  (D. 10) 


Now,  because 


Ep  ‘  OptQl  -  Cp  [R] 


(D.  11) 


and  using  Equations  (D.6)  and  (D.7),  we  have 


=  J 


y2=yl  y2=yl 


Pi 


(D. 12) 


Because,  by  definition  (see  Equation  A. 8) 


C  =  J  . 
P  Pi 


(D. 13) 
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then  E  =  C.  =  C  when  =  p,  .  Also 
P  P  P  2 


H  =  G  +  J1? 

p  p  p2 


(D. 14) 


and  because  equals  zero  whan  J^C^)  =  then 


H  I  *  G  - 

p  1  p  1  pl 

VP1  y2=yi 


(D.15) 


Therefore 


:  *  E  =  G  =  H  =  J  \ 

p  ’  P  P  P  1  p  3 

U2=Pi  y2=VJl  y2=Ul 


(n *  lb) 


We  have 


fp  -  «p  -  %  i*i 


(0.17) 


and,  using  Equation  (D,16) 


FP  I  - 0 

M2=yl 


(D. 18) 


can  be  found  using  the  following  equation 


qJ<53)  qJ<{3) 

F  »  E  — —  -  G  -f — —  +  D 

P  "  P  PJ«3)  P 


(D. 19) 


and  by  using  Equations  (D.16)  and  (D.17) 


Vi 


( n. 20) 


In  addition 


B  *  J  :  +  D  IS]  -  G  [Z] 
p  pl  P  P 


Using 


1 imit  ( Z ] 

U2=U] 


%«2> 

rJ<V 


limit  [S]  =  1 
U2=U1 


,, 

Pl 


and  Equations  (D.16)  and  (D.20),  is  is  easily  shown  that 


B 

P 


0 


Thus ,  we  have 


n  =  B  -  F  =0  when  U9  -  P1 
p  p  P  z 


The  coefficient  Ap  was  expressed  as 


A 

P 


B 

P 


+  C 


P 


(D. 21) 


(D.22) 


(D. 23) 


(D.2A) 


(D . 25) 


(D.26) 


(D. 27) 


56 


A 

P 


=  C 


U2  yi 


P2-Ul  PpKj) 


,X  fei 
pl  eiKi> 


(D. 28) 


Now,  Equations  (30)  reduce  to 


*** '  k 


A  * 

p  ! 

P2= 

=yl- 

pJ(C)  pJ(v> 


(D.29) 


Vtl  A(|jIV  Atv 


y2=Ul 


QX(Q  p\v)  (D.  30) 
P  P 


Comparing  Equations  (D.29)  and  (D.30>to  Equations  (B.la)  and  (B.lb),  which 
are  repeated  for  convenience,  and  using  a  primed  notation 


oo 

E 

T'lr:  1 


Ap'  rj(5>  pJ(v) 


(0.31) 


whe  re 


Fp  QJ(?)  PJ(V) 


P 


F' 


w 
p  pS} 


(0.32) 


(0.33) 


F  ' 
P 


£(v 


pa(f  ) 

p  1 


(0.34) 


k 


V 


Thus 


A  =  A  with  J  (6)  -  J  , (0) 

p  ?  ,  P  P1 

^2  W1 


f;  -  cP  |  with  ye>  =  v0) 

U2=VJ1 


Therefore,  in  the  limit  y0  =  y.  and  with  J  o(0)  -  0, 

Z  i  pZ 


a,  _  =  a;t 
\\>i  tyi 


%II  “  Al|)III  \iv  \v  \ll 


(D. 35) 


(D . 36) 


(D*  37) 


(D. 38) 


;  ******* 
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